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Abstract 

The purpose of this work is to return, with a new observation and rather un- 
conventional point of view, to the study of asymptotically flat solutions of Einstein 
equations. The essential observation is that from a given asymptotically flat space- 
time with a given Bondi shear, one can find (by integrating a partial differential 
equation) a class of asymptotically shear-free (but, in general, twistiing) null geodesic 
congruences. The class is uniquely given up to the arbitrary choice of a complex ana- 
lytic world-line in a four-parameter complex space. Surprisingly this parameter space 
turns out to be the if -space that is associated with the real physical space-time un- 
der consideration. The main development in this work is the demonstration of how 
this complex world-line can be made both unique and also given a physical meaning. 
More specifically by forcing or requiring a certain term in the asymptotic Weyl tensor 
to vanish, the world-line is uniquely determined and becomes (by several arguments) 
identified as the 'complex center-of-mass'. Roughly, its imaginary part becomes iden- 
tified with the intrinsic spin-angular momentum while the real part yields the orbital 
angular momentum. 

One should think of this work as developing a generalization of the properties of 
the algebraically special space-times in the sense that the term that is forced here to 
vanish, is automatically vanishing (among many other terms) for all the algebraically 
special metrics. This is demonstrated in one of several given examples. It was, in 
fact, an understanding of the algebraically special metrics and their associated shear- 
free null congruence that led us to this construction of the asymptotically shear-free 
congruences and the unique complex world-line. 

The Robinson- Trautman metrics and the Kerr and charged Kerr metrics with 
their properties are explicit examples of the construction given here. 
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I. Introduction 



When discussing or studying the general properties of a set of differential equa- 
tions, it is usually best to try to adopt a coordinate system to the symmetries of 
the manifold in question. For example general properties of Maxwell's equations are 
best studied using standard Minkowski coordinates. If however one is looking for 
solutions with specific properties, symmetries for example, it is clearly best to adopt 
coordinates closest to those of the solutions. Spherically symmetric solutions are most 
easily studied in spherical coordinates. In general relativity, when studying general 
properties of asymptotically flat space-time^ the most used coordinate/tetrad sys- 
tem, in the neighborhood of future null infinity, 3 + , are the Bondi coordinate/tetrad 
systenPwhich is associated with the asymptotic symmetry group, i.e., with the Bondi- 
Metzner-Sachs groupP On the other hand, if particular asymptotically flat solutions 
are to be found, special non-Bondi coordinates and/or tetrads are often used. One 
sees this, for example in the study of the Robinson- Trautman metrics or the twisting 
type II metrics where both different coordinates and different tetrad systems are often 
used. 

In this work we will show that for every specific asymptotically flat space-time 
there is a non-Bondi coordinate/tetrad system made-to-order for that solution that 
simplifies the discussion of the solution and brings out its properties more clearly than 
does the Bondi system. Often one can use the Bondi coordinates with the non-Bondi 
tetrad, though for certain situations, discussed later, both the non-Bondi coordinates 
and tetrads are preferred 

To approach this issue, in Sec. II, we give a brief review and survey of the relevant 
properties of 3 + that are needed here. This includes the possible use of two different 
types of coordinates, Bondi and NU coordinates™ and two different types of null 
tetrad systems, (1) those based on a surface forming [twist-free] null vector field near 
'infinity' and (2) those based on a twisting null vector field. In Sec. Ill, we discuss 
what are the different relevant functions that are defined and used on 3 + , (for any 
given space-time), that are used in the different 3 + coordinate/tetrad systems. It 
will be seen that depending on the different choices of coordinate/tetrad system, 
the descriptions of the same physical situation can be very different and that large 
simplifications are available by the appropriate choices. In Sec. IV we will display the 
Bianchi Identities for the different functions in the different coordinate and tetrad 
systems; all quite different from each other. In Sec.V we discuss how one tries to 
extract physical quantities from the variables on J + . 

The major development reported in this work is given in Sec. VI. We show that for 
any asymptotically flat space-time there is a coordinate/tetrad (non-Bondi) system 
that can be adapted to the specific space-time that allows a considerable simplifi- 
cation in the relevant equations. Furthermore it allows us to give a definition of 
angular momentum (orbital and spin) that is geometric and hence independent of the 
choice of coordinates and thus is supertranslation invariant. The point of view that is 
adopted is as follows: an arbitrary asymptotically flat space-time is first completely 
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described in a Bondi coordinate/tetrad frame by the freely chosen characteristic 
data, the asymptotic shear, a, and some (irrelevant for us) initial data. { We have 
tacitly assumed throughout that the space-time can be described or approximated by 
real analytic functions.} Keeping the Bondi coordinates, we now change the tetrad, in 
a specific way: we go from the Bondi tetrad which has non-vanishing asymptotic shear 
(our free data a) to a new tetrad that has (in general) twist but is asymptotically 
shear free. All the information that was in the shear is thus put into the geometry of 
the twisting shear free tetrad. The relevant function to describe the twisting tetrad 
is a complex (stereographic) angle at each point of J + that is denoted by £p + ]. This 
process is reversible. Given L, the a can be reconstructed. The new feature is that 
the shear-free twisting congruence, given by L, is not uniquely determined by the 
original shear. The remaining freedom in L is given by four complex functions of a 
single complex variable which can be thought of, locally, as a (complex) curve, £ a ( T )> 
in C 4 . [Specifically and more surprising, it turns out to be a curve in H -SpaceP^ 
At this point it is a completely arbitrary curve, i.e., any choice of the curve is related 
to the same shear. We show that there is a natural condition that can be imposed 
on the curve that makes it unique (up to initial data) that leads us to define the 
spin angular momentum and the center of mass from the real and imaginary parts 
of the complex curve. These definitions are geometrically invariant under Bondi su- 
pertranslation and reduce to the standard definitions for the Kerr metric. When the 
asymptotic shear is pure 'electric type' the spin vanishes, i.e., the spin arises from the 
'magnetic' part of the shear. Unfortunately to do this construction exactly is quite 
difficult and one must resort to approximations to obtain approximate explicit results. 
Fortunately there is a straightforward approximation method which is summarized in 
an appendix. 

In Sec. VII we give and discuss several examples. They include the Robinson- 
Trautman metrics, the type II twisting metrics and space-times obtained from the 
class of shears referred to as the Tod-Sparling^ shears. 

Though everything that is said and done here could have been stated in the 
language of the conformal compactification of space-time, we will retain the physical 
space-time description. From this point of view, the null boundary is reached by 
letting our null geodesies go to future infinity, i.e., by letting the geodesic affine 
parameters, r, go to infinity. 

We assume that the reader is familiar with spin-coefficient formalism and the 
associated spin-s harmonics though for completeness we have included in Appendix 
C several of the useful relations. 

In addition for comparison with other works, we mention the awkward convention 
adopted here for historic reasons^ and from frequent usage, that our radial coordinate 
r differs from the more conventional choice of r c , i.e., by r = \/2r c and our ub = 
u c / \f2. This leads, in the Schwarzschild case, to the unconventional form 
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II. Review of Scri + : the Future Null Boundary of Space- 
Time 

A. Scri coordinates 

Future null infinity, often referred to as 3 + , is roughly speaking the set of end- 
points of all future directed null geodesies.^ It is usually given the structure of S 2 xR, 
a line bundle over the sphere. This leads naturally to the global stereographic coor- 
dinates, (C, Ci ); f° r the S 2 and the u labeling the cross sections or global slices of the 
bundle. From the point of view of the space-time, 3 + is a null surface and the null 
generators are given by ((, (, ) = constant. 

There is a canonical slicing of 3 + , ub = constant, referred to as the Bondi slicing, 
that it is defined from an asymptotic symmetry inherited from the interior, called the 
Bondi-Metzner-Sachs group. Any one Bondi slicing is related to any other by the 
supertranslation freedom; 

u B = u B + a(C,0- (2) 

with a(C,C) an arbitrary regular function on S 2 . 

All other arbitrary slicings, (u = t), often called NU slicings, are given by 

u B = G(r, CO- (3) 



B. Null Tetrads on Scri 

In addition to the choice of coordinate systems there is the freedom to chose the 
asymptotic null tetrad system, (l a ,m a ,m a ,n a ). Since 3 + is a fixed null surface with 
its own generators, with null tangent vectors, [say n a ], it is natural to keep it fixed. 
The remaining tetrad freedom lies, essentially, in the choice of the other null leg, l a . 
Most often l a is chosen to be orthogonal to the Bondi, ub = constant, slices. In this 
case we will refer to a Bondi coordinate/tetrad system. Relative to such a system 
any other tetrad set, (l* a , m* a , n* a ), is given by a null rotatiorP about n a , i.e., 

j*a = i* + bw? + brn a + bbn a , (4) 

m* a = m a + bn a , 

n* a = n\ 

b = -L/r + 0(r- 2 ). 

The complex function L(ub,C->C->), t ne (stereographic) angle between the null 
vectors, l* a and l a , is at this moment completely arbitrary but later will be dynamically 
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determined. Depending on how L(u B ,C, Cj ) i s chosen l* a might or might not be 
surface forming. We will abuse the language/notation and refer to the tetrad systems 
associated with l* a as the twisting-type tetrads as distinct from a Bondi tetrad even 
when l* a is surface forming. 

The function L(ub,(,C-> ) an d its choice will later play the pivotal role in this 
work. It will be seen that the vector field l* a can be constructed by an appropriately 
chosen L(u, C, C> ) so that it is asymptotically shear-free. 

III. Quantities Defined on 3 + for Any Given Interior Space- 
Time 

Using the spin-coefficient notation for the five complex components of the asymp- 
totic Weyl tensor, we have, from the peeling theorem,^ that 
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where the five, ifj 4 , i(j 2 , ip^ if; , are functions defined on 3 + . For a given space-time 
their explicit expressions depend on the choices of both coordinates and tetrads. 
Their evolution is determined by the asymptotic Bianchi identities and the type of 
characteristic data that is given. 

We adopt the following notation: expression written in the Bondi coordinate/tetrad 
system will appear without a star,'*', e.g., V-L while for Bondi coordinates with a 
twisting tetrad a '*' will be used, e.g., ip^ . In the case of NU coordinates and tetrad 
a double star will be used, e.g. ^2* < \ e t c - 

Using Bondi coordinates and tetrad, the free characteristic data is given only by 
the (complex) asymptotic shear, 

a = o-(u B ,(,(), 

while in Bondi coordinates, but with a twisting-type tetrad the free functions are 

a* = a*(u B ,(,0LL(u B ,(,0 

with a*(u B ,(, () and L(u B ,(, C) carrying the same (redundant) information as did 
the Bondi a{u B) (, (). 
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In the case of NU coordinates and tetrad, the free data is given by 

a** = a**(T,CX)SzV(T,CO 

with 

V(rXX)=du B /dr = d T G(rXX)- 

The important point is that in both cases, the Bondi coordinates with twisting 
tetrad and the NU coordinates and tetrad, all the information that was in a(u, £, Q is 
transferred to the new variables, (<r*& L) or (<r**& V). Later we will show that all the 
information can be shifted into an appropriately chosen L(ub,CiC) with a vanishing 
a* and that in certain special cases (a pure 'electric' type a) all the information can 
be shifted into the V(r, (, () with vanishing a**. 

The relationship between the (^°, ify, ipl, i/jq) given in a Bondi tetrad and the 
(ijj'f, ^3°, V , 2°; i ^0°) °f a twisting tetrad is 

V>o° + 4LVi° + 6L2 ^2° + 4 ^ 3 ^3° + L ^t (5) 
+ SLipf + 3L 2 t/;* + L 3 ^* (6) 

r 2 ° + 2Lr 3 ° + L 2 i,f (7) 
r 3 ° + Lipi° (8) 

<■ (9) 

This is used, in the next section, going from the Bondi version of the Bianchi 
identities to the twisting version. To go to the conventional NU tetrad a further 
scaling is neededP 



< = 
= 

V>2 = 
V>S = 
V>2 = 



IV. The Asymptotic Bianchi Identities 

In each of the different coordinate/tetrad cases there are relations and differential 
equations relating the functions defined on 3 + . Specifically, the ^3 & -0°, in each of 
the cases, are given explicitly in terms of the functions, a, L or V. The dynamics lies 
in the differential equations for the remaining ip's, .i.e., (V'o > > ^2 ) ■ 

Explicitly, in a Bondi coordinate and tetrad system, they are: 



Vg- = -^i + 3a^ (10) 

< = -^° + 2^ (11) 

^0- = _^o + ^o (12) 

ij° 3 = W (13) 

< = -<r (14) 

1^2 — 1P2 — f 2 o" — d 2 a + aa' — aa' , (15) 
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y = ^\ + Q 2 W+aW (16) 
= (17) 

For later use we display the same set of equations in both a Bondi-coordinate/twisting 
Tetrad system and in NU-coordinate/tetrad system but where the simplification 

a* = and a** = 

was used. The justification for this simplification is given in the next section. 

Remark 1 For a Bondi Coordinate/Twisting Tetrad system, with a* = 0, we have, 
from Eq. |3J) ; that 

(18) 
(19) 
(20) 
(21) 
(22) 

(23) 

(24) 
(25) 

Remark 2 Eqs. Ho}) and $%\) are identical even though they appear different. This 
can be seen by using Eq. (H). 

For a NU coordinate/tetrad system, with a** = 0, after rescaling the if)^s with 
appropriate powers of V(r, (, £), 

= -^**° + 3yVC° (26) 
^** 0/ = _^*° + 3yVr° (27) 

^* 0/ = -^* 3 *° + 3—^*° (28) 
V»3*° = 5*3 logV (29) 
^l* = "^y (30) 

^r°-?r° = (31) 

In the first two cases, the dot indicates u^-derivatives, while in the NU case the 
prime means the r derivative. 

We want to emphasize that we have not lost any generality when going to the 
Bondi/Twisting tetrad but the NU equations are valid only when the original Bondi 
shear o~(ub, C> C) was pure 'electric'. In other words all the information in an 'electric' 
shear can be put into the real V(ub, C> 0- Details will be given in Sec. VI. 
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V. The Physical Quantities 



Most often one extracts or attempts to extract the physical quantities by trying 
to identify them in the midst of this set of functions at 3 + by integrals over the ^'s and 
complicated combinations of the a. In general this a difficult and frequently ambiguous 
task - largely because of the Bondi super-translation freedom. Only in the case of the 
identification of four-momentum has it been done successfully and unambiguously. In 
Bondi coordinates and tetrad, one defines the mass aspect $?(v,b,C,C), which is real 
from Eqs.(HH) and (fT7|). by 

= + xPa + aa. (32) 

The four-momentum is extracted from the I = & 1 coefficients of the spherical 
harmonics. When Eq.(|12j) is expressed in terms of $f it takes the very simple form 

= a a (33) 

which, when integrated over the sphere, becomes the Bondi mass loss equation. These 
results can be transformed and reexpressed via the twisting tetrads or in the NU 
coordinate/tetrad system. Explicitly, Eq. (j2Ti|) becomes identical to Eq. fJHH]) except 
that now 

Unfortunately there is no agreed upon asymptotic definition of orbital or spin 
angular momentum. 

In the following section we will address this issue from a totally different point 
of view. Up to now the physical quantities were assumed to be contained in the 
asymptotic Weyl tensor and shear and constructed by a kinematic process. 

We propose that rather than the kinematic process, there is a dynamic process 
involving the known functions on 3 + and only after it is solved can the orbital or spin 
angular momentum be found. 

VI. A New Idea - Dynamics on 3 + 

We first will show that there are special choices of L(ub, (, C) leading to the twist- 
ing tetrads for which the shear, a* = 0. To find this family of functions L(ub,(X)i 
a differential equation must be solved. The solution, however, is not unique. The 
freedom in the solution is four arbitrary complex functions of a single complex vari- 
able, i.e., it is an arbitrary complex curve in a four- dimensional parameter space. 
[The parameter space is the well-studied if -space.] Note that all the information in 
the Bondi characteristic data, ct(-Ub,C,£), will have been shifted to the L(ub,(X)- 
Given an L(ub,(X) with any one of these curves, we could go backwards to re- 
cover the original Bondi o~(ub,(i()- The dynamics lies in the unique determination 
of this curve from other considerations which are describe later. The real part of 
this curve determines the center of mass (or orbital momentum) while the imaginary 
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part yields the spin-angular momentum! 9 * 10 * 1 ^ These quantities are invariant under 
supertranslations. 

We begin by observing that if we start with a Bondi coordinate/tetrad system 
with a given shear a(u B , C 0? then [after an unpleasant calculation^ the shear of the 
new null vector l* a after the null rotation, Eqs.(HJ), is related to the old one by 

a*(u B X,0 = v(uB,(,0- *L-LL\ 

The function L(u, £, £) is then chosen so that the new shear vanishes, i.e., it must 
satisfy 

dL + LL = a(u B ,(,0- (34) 

Remark 3 The special case, dL + LL' = 0, played an important early role in the 
development of twistor theory, leading immediately to the relationship in flat space 
between shear-free null geodesic congruences and twistor theory. It leads us to conjec- 
ture that Eq. \3$ might play a role in some asymptotic form of twistor theory. 

Though Eq.fjHH) is non- linear with an arbitrary right-side and appears quite 
formidable, considerable understanding of it can be found by the following proce- 
dure: 

Assume the existence of a complex function 

t = T(u B , CO (35) 
that is invertible in the sense that it can, in principle, be written as 

u b = X(t,(X). (36) 

Then writing 

(37) 



T ' 

and using the implicit derivatives of Eq. (|3^|) , 

1 = X, T T' (38) 
= d {r) X + X' 9T, 

we obtain 

L = 9 (r) X (39) 

Prime means the r derivative and D( r ) means with r held constant. Finally, with 
this implicit view of X and T, Eq . (J34|) becomes 

d 2 (T) X = a(X,C,0- (40a) 
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Remark 4 Note that this equation, derived here from a different perspective, is iden- 
tical to the so-called good cut equationW^ In general it has a four complex dimensional 
solution space that has been dubbed H-space. H-space possesses a complex metric 
which is Ricci-flat with a self-dual Weyl tensor. It formed the basis of Penrose's 
asymptotic twistor space® In its original form it had its origin in the search for com- 
plex null surfaces that were asymptotically shear-free. Here we are looking for real 
null geodesic congruences that are, in general, not surface forming, i.e., have twist, 
but are asymptotically shear-free^ 

Remark 5 Note that using the reparametrization 

r ^ r * = F(r)=T*(n B ,C,C) (41) 

leaves Eq. \Sl\j and the entire construction invariant. This fact will be used later 
for certain simplifications. 

Remark 6 We mention that if the shear in Eq. is pure electric [i. e., cr(u B , C> C) 

= d 2 S(u B , Cj C) with S(u B , (, () a real function] then the associated H-space is flat and 
has a real four- dimensional subspace that can be identified with Minkowski space. We 
return to this issue later. 



As we just mentioned, the solutions to Eq. (|40aJ) depend on four complex param- 
eters, say z a , and can be summarized by 

u B = X(z a ,(,(). (42) 
However, since we are interested in solutions of the form, Eq. (}36|) . i.e., 

u B =X(r,C,C), (43) 

all we must do is chose in, if -space, an arbitrary complex world-line, z a = £ a (r), and 
substitute it into Eq.(JHJ) to obtain the form (JjSj) . Furthermore, it can be seen that 
the solution, (14*31). can be written as 



u B = X(r, C, C) = C(t)UC, + X m) (r, (, fl 
with Xn^2) containing spherical harmonics, I ^ 2 and 

Tu7^ A, c + c . c-c i-cc > ,,,, 
««) = ^( 1 --TTc7' ! TTc?'TTc? ) (44) 

which has only the I = & 1, spherical harmonics. This result follows from the fact 
that in the good-cut equation annihilates the / = & 1 terms. The £ q (t) does 

get fed back into the higher harmonics via the cr(X, £, £). 

Our solution to Eq. ()34|) can now be expressed implicitly by 



L(U B ,C,0 = * { T)X = e(T)rh a ((,0+ f(r)^2)(T,C,C) (45) 
u B = X = e(r)U(,0+X m) (r,CX). (46) 
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Remark 7 We emphasize that the complex world-line, £ a (V), is not in physical space 
but is in the parameter space, H-space. At this point we are not suggesting that there 
is anything profound about this observation. It is simply there and whatever meaning 
it might have is obscure. We however remark that this observation can be extended 
to the Einstein-Maxwell equations where in general the Maxwell field will have its 
own complex world-line. If the two independent complex world lines coincide, [from 
preliminary results and special cases] one obtains the Dirac value of the gyromagnetic 
ratio^ 

Remark 8 To invert the later equation, fiJW , to obtain r = T(ub, C? C)> general, 
virtually impossible. There however is a relatively easy method, by iteration, to get 
approximate inversions to any accuracy. {See Appendix A} 

Remark 9 We have tacitly assumed that the world line is complex analytic in the 
complex parameter r. 

If we expand £°(t) in a Taylor series and regroup the terms with real coefficients 
and those with imaginary coefficients separately, we can write 

r(r)=&(r)+i£(T). 
This decomposition becomes important later. 

Remark 10 The asymptotic twist, £), {a measure of how far the vector field 

l*a ^ s f rom being surface forming} is defined by 

2%L = DT + LL - dL - LL. 

If L is obtained from a shear that is pure 'electric', then the £ a (r) can be chosen so 
that the asymptotic twist vanishes. 

It is from the fact that we can chose L(u,(,() as a solution of Eq. (}34|) . leading 

to 

a* = a** = 0, that is the justification for the form of Eqs. (fT%|) - (f25|) . 

A. Determining the Complex Curve 

The argument and method for the unique determination of the complex curve 
(and eventually the definitions of spin and center of mass motion) are not along 
conventional lines of thought. To try to clarify it, it is worthwhile to take a brief 
detour. 

It is well-known that for a static charge distribution, the electric dipole moment 
is given by the total charge times the center of charge position given in the 'static' 
Lorentz frame. In the case of a dynamic charge it is more difficult since the center 
of charge depends on the choice of Lorentz frame. Nevertheless by having, in any 
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one frame, the center of charge at the coordinate origin, the electric dipole moment 
will vanish. In the case of a single charged particle moving on an arbitrary real 
world line, (the Lienard-Wiechert Maxwell field) the electric dipole moment vanishes 
when the coordinate origin follows the particle's motion or equivalently, the dipole is 
given by the particles displacement from the coordinate origin times the charge. The 
asymptotically defined dipole moment also vanishes when it is calculated (extracted) 
from the I = 1 part of the coefficient of the r~ 3 term of <ft = F a bl a m b , where l a 
is a tangent vector of the light-cones that are attached to the particles world-line. 
Analogously, it was showrPl that the magnetic dipole moment of a charged particle 
can be 'viewed' as arising from a charge moving in complex Minkowski space along 
a complex world-line. It is given by the charge times the imaginary displacement. 
{We emphasize again that the complex world-lines are a bookkeeping device and no 
claim is made that particles are really moving in complex space. The real effect of 
the complex world-line picture is to create a twisting real null congruence.} If the 
light-cone from the complex world-line is followed to 3 + , (equivalent to introducing 
an appropriate asymptotic twist for the null vector l a that is used in </> = F a bl a m b ) 
then both the asymptotic electric and magnetic dipoles vanish. This argument is 
given in considerable detail in the paper.^E^ 

The idea is to generalize this construction to GR. In linearized GR and in the 
exact case of the Kerr and charged Kerr metrics, the mass-dipole moment and the 
spin sit in the real and imaginary parts of the I = 1 harmonic of the coefficient of 
the r -4 part of the Weyl component if)® in a Bondi coordinate/tetrad system. This 
observation is extended to all asymptotically flat vacuum solutions by going to an 
asymptotically twisting tetrad, with an appropriately chosen if-space complex curve 
[see below], such that the new Weyl component, ip] , has a vanishing I = 1 harmonic 
in the coefficient of its r -4 part. We refer to this curve as the intrinsic complex center 
of charge world-line and, rough speaking, the center of mass and the spin angular 
momentum will be identified by its real and imaginary parts. Effectively, by chosing 
the world-line so that the I = 1 harmonic vanishes we have shifted the 'origin' to the 
complex world-line. 

To explicitly apply this idea to our discussion of twisting null congruences, we 
go to the set of equations, Eq. (fTT?|) and (|2*U|). namely 



and require that, in (}4"Tj) . that the I = 1 part of tpf should vanish. This leads (via 
approximations) to a differential equation for the complex center-of-mass world-line. 
It is driven by the original Bondi shear. To see this in detail, first we recall that (J48j) 
can be rewritten as 



-ty?-Li/>?--2i4>« 



(47) 




(49) 
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a = c)L- + i(L 2 r (50) 
$f = m - 2L - L 2 W - D 2 o - oW. (51) 

If we use the result [see Appendix B] that the Bondi four- momentum, P a , can be 
extracted from \I> by 

c 2 



P a (u B ) = (M, -n = -— I m a dS (52) 



we can invert and see that 



*(u B , C, C) = X - X% + ..- = -2- 2 {V2M - 3P%) + * {>2 ... (53) 

c z 



with 



S(C,C) = ^(i,a(C,0), (54) 



c b = l b - % = V2(0, a(C, 0) (55) 

C a = yc , (56) 

Cj a unit radial vector and \I//> 2 containing only harmonics I > 2. 

The time evolution of P a {u) is found by integration Eq. (|49p . i.e., from \1/ = 
J o'a'du or 

P a W = - / aW? a dS. (57) 



If ?/>2°; from Eq. (|51J). is substituted into f)47l ) we obtain after using Eq. (|49jl and 
regrouping, a truly ugly equation that we write as 

- 3 a+ W + 3m = (NL) (58) 
with the non-linear terms, (NL), expressed by 

(NL) = D[Wa] + la W + 3L d 2 W + L 3 a"' + 3LaW (59) 
+3L 2 dW + 3L[2L W + L 2 a" + d 2 a + era']. 

The grouping has the non-linear higher order terms put into (NL) while the 
'controllable' (lower order) terms are on the left side. It is from this equation and 
(|H7j) we find the equations of motion. We first require, as we said earlier, that tp{° 
(and of course ipf' ) have no I — 1 harmonics. From the fact that t) 3 cf, from its spin 
s = 2 structure, has no I = 1 terms we see that Eq. ()58j) reduces to the condition that 

(W + 3L-*) l=1 =(NL)i =1 . (60) 
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It is this equation that we must analyze, or more accurately approximate, by looking 
only at the leading spherical harmonic terms. [See the following subsection for an 
alternative derivation of Eq.(|6L)|h] 

From Eqs. (|53l ) and ([45)1 . considering only the lowest order I = 1 terms we have 

*(UB,C,C) = -2%{^2M-W l Cl ) 
L(u B ,U) = C(r)m a ((X) 
from which it follows that 



G 

9^ = 12— P a m a (61) 
c 2 

3LV = -Q^-C'(r)m a {V2M -3P l a} (62) 

cr 

= -V26^MC-(r)m a + 18^C' P\r)m a c b (63) 
c c 

= -V2Q^MC'(r)m a + 9i^C P b (r)e abcd m c t d + Y{1 ;> 2} (64) 

where we have used the Clebsch-Gordon expansion [see Appendix C] of the product 
m a c b . From Eq.()6U|). we then have 

12^P a m a - V26^MC\r)m a + 9i^ a 'P b (r)e abcdV cf m f t d = (NL) l=1 = GN a m a 



or, using the rescaled « B to the conventional Bondi u c coordinate ^ by ub 

P e = Mf'(r) - i^e'P\r)t abcd r 1 ce t d + D e . (65) 

We have the three-momentum expressed in terms of the mass M and the complex 
velocity, £ a '. Remembering, from Eq. (j25j) that ^ is real we see, from Eq.(|53p. that M 
and P c are real. By defining the real and imaginary parts of £ a ', 

r = (66) 

and taking the real and imaginary parts of ()65|) . (ignoring in this approximation) 

P e = Me R {r) + ^CiP\r)e ab c d V ce t d (67) 
MQ\r) = \e R P\r)e abcd ^t d . (68) 

If we identify the spin-angular momentum (from the analogy with the Kerr met- 
ric) by 

S b = Mfi 
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and £ R = v with a real velocity vector and linearize, we have 



P c = Mv c 
g- = 

M 



which, with 



P a (u B ) = (Mv c y = - J atrTadS, 



yields the equations of motion for the center of mass and spin-angular momentum 
- the equations of motion of a spinning particle with mass and momentum loss via 
radiation — all derived from the Bianchi identities and driven by the original Bondi 
shear. 

Remark 11 Note that from the linearization of Eqs. \61\) and \6§). the coupling be- 
tween the spin and the velocity vector has disappeared. The details of the coupling will 
be investigated in a future note. 

There is an important point to be made here. We know that the complex curve 
we wish to determine is given by the expression 

with the complex paramter r, while in the above expressions for the determination of 
the curve we have treated £ 6 to be a complex funtion of the real parameter u. This is 
justified by realizing that we have been making severe approximations, throwing out 
higher non-linear terms, not using more and higher Clebsch- Gordon expansions and 
from the linearization of the inversion of r = T(u,(,(), Appendix A. In the lowest 
order of the inversion we do have r = u. 

The basic idea we are espousing is that there is, in if-Space, a unique complex 
curve that can be called the intrinsic complex center of mass world-line. It is deter- 
mined from a geometric structure on the physical space-time, the null direction field, 
independent of the choice of coordinates on Scri. In order to obtain - as an approx- 
imation - what one would normally call the orbital and spin angular momentum in 
linear theory one must use Eq.(JHJ and extract by an integration, from the Bondi 
its I — 1 coefficient. The values of this 'normally' defined angular momentum will, in 
general, depended on the choice of the scri coordinates or the Bondi frame. 

Remark 12 If we linearize (decouple) Eg. \65\) we have 

P e = MC\u) 

so that if we treat M and P e as very slowly varying we can integrate it as 
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or, since M and P e are real, 

pe 

= const = initial conditions, 
thus obtaining the usual motion for the center of mass. 

B. An Alternative Route to the Complex Curve 

There is an alternative method of finding the curve that in some sense is more 
direct than the one just given. 

We begin by inverting Eqs.© " © an d then writing 

< = V? - 3i^2 + 3£ 2 ^ - L^l 

and its w^-derivative 

V5° ' = - 3L$j + 3LV§ - L^l) ■ (69) 

We then impose our basic requirement for the determination of the curve, namely 
that the lowest harmonic, i.e., the I = 1 term in and should vanish. Eq.()69jl 
can then be written as 

{Vi'}i=i = {(3^2 - 3L 2 ^3 + ^VS) }i=i + [Tevms]^. (70) 

Eq. (|7()jl . using the Bianchi identities for ip®' and ijj^' and expressions ip 3 and i.e., 
Eqs.(|ll|) - (|14|) . again leads to the complex curve by the following argument. 
Using 

$ =^0 + 5 2- +0 .-. 

instead of ^ an d the evolution equation for Eq. ()33|) . (really a Bianchi identity) 

r = a ct (72) 

Eq.(|70J) can be rewritten as 

{3L * - ^?'}i=i = + [ Terms]^ 2 , (73) 

while the Bianchi identity for ip® ' becomes 

ipi' = — + NL + [Terms] (74) 

In both cases, for our approximation, we have simply grouped the higher har- 
monic terms and the nasty non-linear terms into [Terms]/j>2 and NL. 
By substituting the -0°' from Eq.([74|) into (J75J) we obtain 

(X>* + 3L-*) l=1 = (NL) l=1 , (75) 
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which is identical to Eq.()60|). This relationship, with Eqs. (j45|) and ()46|) . 



L(u B ,(,0 = r(r)m a (C,C)+ f W X ( ^ 2) (r,C,C) (76) 

UB = e(T)U(,()+X il>2) (T,(,0 (77) 

leads, via the argument in the previous section, to the equations for the complex 
curve. 



VII. Examples 

A. The Robinson- Trautman metrics 

We begin with the Bianchi Identities Eqs. (j26j) - (J31|) . which are essentially the 
same as Eqs. (J18j) - (}2l)j) . differing only by the coordinate transformation r =^> ub, given 
by ub = X(t,(,C)j from one set to the other. The tetrad system, as well as the 
stereographic angle L, are the same with both sets. In other words, geometrically, they 
represent the same situation. The L or the l* a describe a surface forming congruence. 
We imposed, in the last section, the condition that the I = 1 harmonics of the tetrad 
component, tp®, vanish. In this example we simply strengthen this condition by 
requiring that 

ifi = 0. 

This condition, from Eq. pTj) forces ip^ to be a function only of r, i.e., ^\ = ^(t) = 
2 v / 2GM(r)/c 2 , so that Eq.tpKj) becomes 

V' 

^* 0/ = _00c)c) log V + 3y^j* 2 *°. (78) 

(The strange numerical factor 2\^2 arises from our use of the radial coordinate r that 
differs from the conventional one by the factory^-) 

By the reparametrization freedom, Eq. (j41|) . a change in the r coordinate can be 
made, i.e., r = F(r*), so that the -02*° i s constant, thus yielding 

— = 0009 log V (79) 

the well-known Robinson- Trautman equation^ describing an asymptotic Robinson- 
Trautman space-time. If we further assumed that ip = this would lead to the 
Robinson- Trautman metric. 
Note that since 

u = x = a^K(c,c)+^wr,c,c) 

v = x' = e'(r%(c,0 + x{ l>2) (T,c,0 
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we see that Eq.(|79j). is a differential equation for V and more specifically the coeffi- 
cients of the I = & 1 harmonics have the form of Newton's second law of motion 

MC" = F a 

so that the Robinson- Trautman space-time describes a gravitational rocket ship. 

If we do have a solution to the Robinson- Trautman equation, i.e., V = V(r, (, (), 
the associated Bondi shear can be easily constructed from the implicit relations 

u B = X(r,C,C) = J V(T,(,()dT (80) 
a(u B ,(,0 = f? T) X(T,C,C). 



B. Twisting Type II Metrics 

We can describe a more general situation by going to Eqs.(JTSJ)-(|13J) and again 
require that the Weyl components 

= 

< = o. 

This leads to the equations 



(81) 
(82) 



with 



tP* 3 ° = W + L&-, (83) 

$ = + 2L W + L 2 W + d 2 a + aa, 
^ = 0, 

that are precisely the equations for the twisting algebraically special type II metrics. 
A detailed analysis of these equations and their generalization to the type II Einstein- 
Maxwell equations is in preparation. 



C. The Sparling- Tod Shear 

George Sparling and Paul Tod 6 found a class of Bondi shears such that the 
good-cut equation ()40a|) 

d 2 X = a(X,(,0 (84) 
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can be integrated exactly. An example of their class is 

, , 7 , (qgm a ) 2 

<T(V>,tX) = 3 — (85) 

u 

with q a the constant vector 

g a = (0,ft). (86) 
The solution to (|H1|) can be given by 

v 2 = X\z a XA) = Z 2 + S 2 (87) 

with 

Z = z a l a (CX), (88) 

s = s a i a ((,c). 

The four z a = (t, x, y, z) are the if-Space coordinates (four constants of integra- 
tion) and the s a (z b ) are four functions of z a determined by the algebraic equation 

s a z b (l a m b - l b m a ) = q a m a . (89) 
Using the identity (see Appendix C for the notation and conventions) 

(T a rh b - mj b ) = ^{t a m b - m a t b ) - ^ie abcd m c t d (90) 

and the ansatz 



s a = S P + Si q a + i S 2 V2e a bc q b z c (91) 
V2e a bc = -e a bcd t d (92) 

one find that the (So, Si, S2) are uniquely determined so that 

s a = — e {z b q b t a - z h t b q a + i e a bcd t d q b z c }. (93) 

z e z e 

Note that though s a , in Eq. ()89j) was defined up to a term proportional to z a , the 
ansatz determined it uniquely 

Eqs.(jnHJ), (jHTj) and (??) determine the solution ub = X(z b ,(X) so when an 
//-Space world-line, z b = £ 6 (t), is chosen we have the null direction field 



£(«B,C>C) 



9 (T) X(r(r),C,C) 

Vz 2 + s 2 . 



ZdZ + SdS 

Vz 2 + s 2 



(94) 
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D. Arbitrary Construction 



For further general examples, we can reverse the process of first giving the shear 
and then finding the L(ub, C> 0- Instead, we can choose an arbitrary, s = 0, function 
with the form 

u B = X(r,C,C) = r(rK(C,C) +X (m (rXX) 

and then define 

L(u B ,CX)= DX = C(T)m a ((,()+ WWr,C,C). 
The associated Bondi shear is then given, implicitly, by 

a(u B ,C,C) = * 2 X {1>2) (t,(,() 

u = C(r)Uc,0+X m) (r, CO- 
VIII. Conclusion 

In this work we have shown that for all asymptotically flat space-times there is 
a hidden structure that must be extracted dynamically from the known asymptotic 
Weyl tensor and its related characteristic data, the Bondi shear, a. This hidden struc- 
ture is a specific field of asymptotically shear-free null directions - or an asymptotically 
shear-free null geodesic congruence. Within the information for the description of this 
direction field is a complex world-line that is defined in the if-space associated with 
the given asymptotically flat space-time. One can try to give meaning to this world- 
line by defining it to be the complex center of mass of the interior gravitating system. 
Though there is no proof or even overwhelming evidence that this assignment is com- 
pletely reasonable, we can ask the question: what physical justification can be given 
for this assignment. On the surface it certainly is strange - where does a complex 
world-line in if-space enter in any direct physical observation. 

We would like to try to give a justification by three different types of argument. 

1. Examples: In several known cases, e.g., the Kerr-metric, the charged-Kerr 
metric, the Lienard-Wiechert-Maxwell fields, the Kerr-Maxwell field, all have a center 
of mass (or center of charge in the Maxwell case) that is identical to the more general 
construction described here. In addition, the construction given here coincides with 
the definitions from linearized theory. 

2. Though it appears to us to be surprising and remarkable, one does get, in a 
first approximation, the conventional equations of motion of a complicated spinning 
gravitating source directly from the asymptotic information - with no details needed 
about the interior source. These equations contain the radiation reaction from the 
gravitational radiation that arises from the Bondi mass loss. 

3. It is possible to ask the question: can one give, at least in principle, an 
observational means of "observing" this complex motion? The answer, we believe, 
is yes - though to really do so is impossible. It involves having a huge number of 
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observers surrounding the gravitating source and looking at all the null rays reaching 
them, then picking out shear-free null direction fields and by angular integration 
finding the complex world-line. This type of argument must be tightened and made 
more precise. 

Work has begun on the application of these ideas to asymptotically flat Einstein- 
Maxwell fields. We remark that it appears (in a preliminary stage) that there are two 
complex world-lines, a complex center of mass obtained from the asymptotic Weyl 
tensor and a complex center of charge obtained form the asymptotic Maxwell field. 
Their imaginary parts give respectively the spin and magnetic dipole moments. If the 
two world-lines coincide then the Dirac value of gyromagnetic ratio follows.^ 

Many years ago, by a rather strange and not at all understood complex coordinate 
transformation applied to the Reissner-Nordstrom Einstein-Maxwell field, the charged 
Kerr Einstein-Maxwell field was first found.^ The work described here appears to 
clarify and even explain this mysterious transformation. 

Nevertheless, there is still much to be understood about the present material, 
i.e., the asymptotic shear-free congruences and the complex H-space curves. Is it 
just an coincidence that it appears to be so connected with physical issues or does 
something deeper lie behind it? 
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XI. Appendix A 

In order to obtain r = T(u B , (, C)? the inversion of 

u B = X(r,C,C) = z a (r)T a (CC) +X m) {r,CC), 

we first note that r can be replaced by any function of r, i.e., the Eq. (J37|) is invariant 
under r = F(r). See Eq.([41j). Thus by taking 

z°(r) = r 

we have 
or 

t = u B - z\T)k(CX) + X m) (r,CX). 
This can be solved by iteration to any order: 

TO = U B 

n = u B - z l (u B )k((,() + x {l>2) {u B ,cX) 

T n = U B - Z*(r n _i)/i(C,C) +X(i22)(T n -l,(,()- 

XII. Appendix B 

Normalization and conventions: Our goal here is to express the function ^ in 
terms of the Bondi mass/momenta and also show how to extract them from the 

We start, in a given frame in Minkowski space, with C a a unit radial space-like 
vector and T a a unit time-like vector given by 

C a = (0, cos sin 6 1 , sin sin ^, cos#) = —C a , (95) 
T a = (1,0,0,0). (96) 

with the null vectors 

l a = — -= (1, — cos0sin0, — sin0sin#, — cos#) = — j=(T a + C a ) (97) 
v 2 V 2 

n a = —= (1, cos sin 6 1 , sin sin 6 1 , cos^) = —=(T a — C a ) (98) 
v 2 v 2 

and 

c a = l a — n a = — \/2(0, cos sin 0, sin sin 0, cos 9) = V2C a 

We let a = (0, 1, 2, 3) = (0, i) and l a = (/ , ^) and = v^2Cj, etc. 
We then define the / = and Z = 1 parts of the mass aspect ^ by 
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Therefore 



* = ip* 2 ° + 2L W + L 2 cT + d 2 a + a? (99) 

* = x-x% + ~- (ioo) 



= Xla-X%la + ... (101) 



^=[ X T a - v^^Wa + xC a - V^CX 



From the integral identities 





4tt 


y = 









C iC jdS = 


3 J 



we have that 



J m a dS = 1= j dS( X T a - VltfCiCa) 



with x 4 = -X 1 



On the other hand, from reference (f^) and a change of notation, " we have that 
the 4-momentum for an asymptotically flat space-time is defined by 



2 (• 
07T(jr 7 

* = -2^2— (M+ (104) 



c 2 



then from Eqs. ([T03|) . ([T02|) and (fTTTTjl . 



p « = (105) 
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We find that: 



X 



-2y/2%M (106) 



c 2 " 



Xl = -6GP t = 6^P l (107) 



c 2 ~ 



and therefore we have 



* = X - X% + .... = -2^(V2M - W%) + .... (108) 



c 2 



XIII. Appendix C 

Relations between the following quantities have been used through this work: 



? = ±±1 -i^- ~ 1 + C h (109) 

2 'l + CC' i + CC' i + cc 



2 i + CC i + CC i + CC 

m = aft = ^(o,I^,^±0,^), 
2 i + CC i + CC i + CC 

ft = v^(l, 0,0,0), 

n a = ft ft- Al C+ l i C ~l 1-C ^ 

2 ' i+_CC' l + CC'i + CC 

c" = T a -n a = y/2(0, -^4, -i-^i, _1 + - C N 

i + CC i + CC i + CC ' 



with the products 



2 

(m a m b -m a m b ) = --e abcd (?t d , 
(rh a Cb - c a m b ) = ie abcd rh c t d . 

m a c& = -(m a c fe - m b c a ) + -(m a c fe + m b c a ) 

= i^ abcd m c t d + {l = 2)Terms 



and 

e = — 

v/2 



^abc /—^abcdt ■ 
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